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Abstract. We present the results of 5 GHz VLBI obser- 
vations of a compact steep spectrum source 3C 138. The 
data are consistent with the western end being the loca- 
tion of the central activity. The observed offset between 
different frequencies in the central region of 3C 138 can 
be accounted for by a frequency dependent shift of the 
synchrotron self-absorbed core. Our new measurements 
confirm the existence of a superluminal motion, but its 
apparent velocity of 3.3c is three times slower than the 
reported one. This value is consistent with the absence of 
parsec-scale counter-jet emission in the inner region, but 
seems still too high to allow the overall counter-jet to be 
seen in terms of Doppler boosting of an intrinsically iden- 
tical jet. Either an interaction of jet with central dense 
medium, or an intrinsically asymmetrical jet must be in- 
voked to reconcile the detected superluminal speed with 
the observed large scale asymmetry in 3C 138. 
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1. Introduction 

The radio source 3C 138 (= 4C 16.12 = J0521-hl638) is 
identified with a quasar with 77i,„=18.84 and z=0.759 (He- 
witt & Burbidge 1989). It is a compact, powerful quasar 



observations (Akujor et al. 1995| and references therein; 



with a well-defined turn-over in the spectrum at about 
100 MHz and a steep high frequency spectrum of 0.65 
(S (X making it a prototype of the com pact steep 

spectrum (CSS) radio sources (Fanti et al. 199C). Its spec- 
trum remains straight up to 22 GHz, supporting t hat n o 
dominant flat spectrum core exists (Kameno et al. 1995|) . 

The radio structure of 3C 138 consists of a core with 
a bright jet and compact lobe to the east and a fainter, 
more diffuse lobe to the west from the VLA and MERLIN 



Liidke et al. 1998). High resolution VLBI observations 
(Geld zahler et al. |1984|; Fanti et al. |1989|; Nan et al. 



1992 ; Dallacasa et al. 1995| ; Cotton et al. 1997) reveal sev- 



eral knots within the main jet emission extending about 
400 mas in a position angle of 65° and a low brightness 
counter-jet about 250 mas away in the opposite direction. 
There are two kn ots d ominating the central 10 mas re- 
gion. Fanti et al. ( |1989|) claimed on the basis of ap parent 
size that the eastern knot is the core. Cotton et al. (1997) 
argued that the western one is the location of the nuclear 
core mainly because of its very weak polarization (< 0.4%) 
compared to a peak polarized intensity of 3.5% for the 
eastern component at 5.0 GHz. Cotton et al. ( 1997 ) also 
suggested a relative proper motion of 0.27 mas yr~^ be- 
tween the two knots, corresponding to an apparent veloc- 
ity of 9.7c (assuming Hq = 65 km s~^ Mpc~^ and qo = 0.5, 
hereafter). 

In general, the radio emission from CSS sources is 
thought to be relatively free of beaming effects, consis- 
tent with the fact that CSS sources are the least variable 
of compact extragalactic radio sources, and CSS sources 
are rarely found to be superluminal (O'Dea 1998 ). When 
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combined with the observed two-sided jet structure, such 
a superluminal motion in 3C 138, if confirmed, will strictly 
confine the jet velocity and its direction to the line of sight 
of the observer. We report on 5.0 GHz VLBI observations 
of 3C 138 as an attempt to confirm its proper motion and 
probe its compact radio core area. 



2. Observations and data reduction 

The VLBI observations were carried out at 5.0 GHz on 
November 8, 1997 with the enlarged European VLBI 
Network composed of Effelsberg (100-m), Hartebeesthoek 
(26-m), JodreU MK2 (25-m), Medicina (32-m), Nolo 
(32-m), Onsala (25-m), Shanghai (25-m), Torun (32-m), 
Urumqi (25-m) and Westerbork (25-m). Sources were ob- 
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Fig. 1. Tapered 5 GHz VLBI image of 3C 138 restored with a 10 mas circular beam. Contour levels are drawn at —3.6, 
3.6x2" mJy beam-^ (n = 0, 1, 6). The peak flux density is 274 mJy beam'^ The RMS noise is 1.2 mJy beam . 



served in a snapshot mode, i.e., 5-12 thirteen-minute scans 
spreading across a wide range in hour angle to obtain a 
good u-v coverage. The left-circular polarized radio signals 
were recorded in Mark III mode C with a total bandwidth 
of 14 MHz (7 IFs and 2 MHz per IF band) at each station. 
The correlation of the data was performed on the MPIfR 
MK HI Correlator in Bonn, Germany with an output av- 
eraging time of 4 s. 

Then, a priori visibility amplitude calibrations for each 
antenna were applied using the antenna gain and the sys- 
tem temperatures measured during the observations. Af- 
ter the data were globally fringe-fitted in AIPS (Schwab 
& Cotton 1983), the subsequent data inspection, imaging 
and modelfitting were carried out with DIFMAP (Shep- 
herd 



3. Results 



3.1. Overall structure 



1997) 



Because of the large separation (^400 mas) between 
the eastern compact lobe and central core region in 
3C 138, special attention was paid to the smearing ef- 
fects of time averaging and bandwidth in order to image 
an undistorted full field. The 2 MHz data in each IF were 
not averaged through 7 IFs, and the averaging time was 
only 12 s to minimize smearing effects. 

Images were produced with the usual iterations of 
cleaning and phase self-calibration. Only a constant gain 
factor for each IF band of each antenna was implemented 
in the later stage of imaging. These gain corrections are 
consistent with those obtained for other sources (Shen 
et al. in preparation). To improve the sensitivity to the low 
surface brightness regions, the visibility data were natu- 
rally weighted and tapered with FWHM = 30MA. 



The overall morphology of 3C 138, as shown in Fig. ^ 
resembles past VLBI images. There are two notable emis- 
sion regions: the central core and the eastern lobe, which 
are separated by 400 mas at a position angle of 70°. In 
between them are several discrete emission components. 

The total cleaned flux density in the image (Fig. ^ is 
~0.85 Jy, which is the same as the correlated flux density 
on the shortest baseline, Westerbork-Jodrell MK2, but far 
less th an th e values of 3.24 Jy a nd 3. 58 Jy by Dallacasa 
et al. ( |1995D and Cotton et al. ( |l997[ ), respectively. The 
single-dish measurements from the University of Michi- 
gan Radio Astronomy Observatory (UMRAO) show a very 
stable total flux density of 3.80 Jy at 4.8 GHz with a varia- 
tion less than 6% over past 10 years. The source is heavily 
resolved on baselines longer than 3 MA (see Fig. 2c in 
Fanti et al. 1989). So, the main reason for the shortage 
of cleaned flux density in our image is the lack of shorter 
baseli nes fi ke MERLIN or VLA-Pie Town (cf. Dallacasa 
et al. 1995). As a result, the more extended, lower surface 
brightness regions seen in previous images (e.g. Fanti et al. 
1989) are not well recovered in our image, which causes 
discontinuity in the main jet emission towards the eastern 
lobe. The connection of separated emissions (including the 
core region), however, exhibits a sign of oscillation seen 
from other observations. No counter-jet was significantly 
detected in our image though there are some very weak 
diffuse emissions at about 3 times noise level around the 
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Table 1. Gaussian model for the central region of 3C 138 





Right Ascension (mas) 

Fig. 2. Full resolution 5 GHz VLBI image of the 
central core region of 3C 138. The resolution is 
1.2 masx 0.92 mas at a position angle of 37° (indi- 
cated at the lower left corner). Contour levels are drawn 
at -3.6, 3.6x2" mJy beam^i (n = 0, 1, 5). The 
peak flux density is 134 mJy bcam^^. The RMS noise 
is 1.2 mJy beam^^. 

count er-jet region of the 1.7 GHz image by Cotton et al. 
(119971) . The RMS noise in the image is 1.2 mJy/beam, 
about twice the expected thermal noise for the total on- 
source time of 1.3 hr. 

Two compact components (labeled F and G after Fanti 



et al. 1989) in the eastern lobe can be clearly seen. They 
are aligned nearly orthogonal to the jet axis, which may 
be interpreted as a working surface of the jet due to its 
interaction with the surrounding medium. 

3.2. Central core region 

For the purpose of identifying the center of activity, a 
full resolution image is displayed in Fig. |[ This image 
was made from the same calibrated data as Fig. |l|, except 
that the baselines shorter than 30 MA were removed and 
the uniform weighting was adopted to ensure a high res- 
olution. To yield a quantitative description, modelfitting 
to both amplitudes and phases in the calibrated visibility 
data was applied. Because of the complicated overall struc- 
ture and the limited observational data, we fit a model to 
core region only while keeping the cleaned components for 



the remaining structure (cf. Finer et al. 200C). In this way, 
we can minimize the effect mainly from the components F 
and G. We also tried a modelfitting to the longer baselines 
data corresponding to Fig. ||. We even made an image of 
the central region only with a 1-minute averaging time 
and then fit the visibility data. All these results are quite 
consistent. In Table 1 is presented the best fitting model 
consisting of three circular Gaussian components. 

The sizes of components C^ and Gw are about the 
same, but Cg is much brighter. The flux density of C^; 



is stable over past 12 years (cf. Akujor et al. 1993; Dalla- 
casa et al. 1995). The two point spectral indices (as ocHz) 
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S: the flux density in Jy; (r, 6): the distance and position an- 
gle of each component with respect to component Cw in mas 
and degrees, respectively; a: the diameter (FWHM) of circular 
Gaussian component in mas; afioGHz- two point spectral index. 
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Fig. 3. The apparent separation of two central compo- 
nents Gw and C_e in 3C 138 as a function of the observing 
epoch. These measurements are made at 5.0 GHz except 
one (denoted by a open triangle) at 8.6 GHz. Symbols 
are explained in the text. Also shown is a fitting line (see 
text). 



obtained by combining the high resolution VLBA results 
at 8.6 GHz (Fey & Chariot 1997) are almost the same for 
both Ge and Gw (see Table 1). Component Ge2 is less 
well constrained because it is much too resolved. Actually, 
C E2 emission has appeared in several previous 5 GHz im- 
ages (e.g. Fig. 13 in Cotton et al. (1997)), but was first 
fitted at 8.6 GHz by Fey & Chariot W99%. 



4. Discussion 



4.1. Proper motion 



There are three published 5 GHz VLBI measurements 
made at epochs 1985.5 (Fanti et al. |1989D , 1 989.7 2 (Dal- 
lacasa et al. |1995D and 1994.97 (Cotton et al. |1997| ). Fanti 
et al. (1989) resolved the central region with two com- 
pact components, but did not present the separation at 
the time. Dallacasa et al. ( |1995D made the first measure- 
ment of the apparent separation of 5.0 mas. Cotton et al. 



(1997) estimated the separation of 3.9 and 6.5 mas for 



epochs 1985.5 and 1994.97 respectively, and obtained an 
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averaged rate of separation of 0.27 mas yr^^ from all three- 
epoch observations (see Table 1 and Fig. 14 of Cotton 
et al. (1997)). A separation of 6.5 mas at epoch 1994.97 
is larger than that of 5.76 mas at a later epoch 1997.86 
in Table 1. This would be indicative of a back flow in 
the jet motion, which is very unlikely if not impossible. 
The 8.6 GHz VLB A observations (Fey & Chariot [1997D 



measured a separation of only 5.4 mas at epoch 1995.78, 
which is also much less than that at epoch 1994.97. To fur- 
ther investigate this, we measured the peak-to-peak sep- 
aration from images directly, and then compared them 
with those modelfitting results. As a result, we got an 
excellent agreement between the modelfitted separation 
and the peak-to-peak separation at epochs 1989.72 and 
1997.86 at 5.0 GHz and at epoch 1995.78 at 8.6 GHz, but 
the discrepancy at epochs 1985.5 and 1994.97 was as large 
as 0.5 and 0.9 mas, respectively. This difference is about 
one-half and one-quarter of the corresponding beam size, 
and quite difficult to understand for us at present. 

Therefore, we decided to estimate the proper motion 
by using two epoch (1989.72 and 1997.86) 5 GHz data 
(denoted by crosses in Fig. ^) only. As shown in Fig. ||, 
this produces a relative proper motion of 0.093 mas yr 
between C^ and Cw, which corresponds to a superlumi- 
nal speed of 3.3c. We used other two 5.0 GHz observations 
to evaluate its uncertainty. The peak-to-peak separation 
(denoted by open squares in Fig. |) of 4.4 and 5.6 mas 
at the corresponding epochs 1985.5 and 1994.97 gave an 
uncertainty of 5%, while the published points (denoted by 
filled squares in Fig. ^ would introduce an uncertainty of 
16%. This proper motion is much slower than the reported 
0.27 mas yr"-'^ (Cotton et al. 1997 ). The ongoing monitor- 
ing observations confirm that the knot is definitely moving 
much slower than suggested by the earlier data (Bill Cot- 
ton, private communication 2000). 



4-2. Synchrotron self- absorption 

There is a relative position offset of 0.17 mas at epoch 
1995.78 for Ce between the position interpolated from 
5 GHz data (fitting line in Fig. ^) and the measure- 
ment made at 8.6 GHz (Fey 
by a open triangle in Fig. H. 



Chariot 1997) denoted 



Such a frequency depen- 
dence of the observed position can be explained by syn- 
chrotron self-absorption (Lobanov 199^ ) . To estimate this 
position offset using eq. (11) of Lobanov ( |1998D , we have 
searched the literature for its radio spectral energy dis- 
tribution. As a result, the total synchrotron luminosity 
of 1.5x10"*^ erg s^^ is estimated. A superluminal speed 
of 3.3c implies a maximum angle 9 to the line of sight, 
a minimum jet speed v (= /3c) and a minimum Lorentz 
factor 7 of 34°, 0.96c and 3.5, respectively. The projected 
opening angle 4>o of 5.8° is measured from Ce component 
size over its distance to C^ at 8.6 GHz. These values pre- 
dict a maximum position offset of 0.22 mas between 5.0 
and 8.6 GHz under the assumption of energy equipartition 



between radiating particles and magnetic field. The upper 
limit is mainly due to the uncertainty in the determina- 
tion of the viewing angle 9 and the resultant lower limits 
to both f3 and 7. Thus, such an observed offset in the cen- 
tral region of 3C 138 can be well explained as a frequency 
dependent shift of the synchrotron self-absorbed core. As- 
suming a conical jet geometry in which the magnetic field 
decreases linearly with its distance to the core, we can 
further derive the magnetic field in the 5 GHz VLBI core 
region of 0.12, 0.02 and 0.01 G for components Cw, Ge 



and Ce2, respectively (cf. Lobanov 1998) 



Recently, a significant position offset has been reported 
for another superluminal CSS source 3C 216 (see Fig. 3 of 
Paragi et al. 2000). These may suggest that the core re- 



gion in CSS sources could also be described by ultracom- 
pact jet models, which have been successful in interpreting 
the central region of radio-loud active galactic nuclei (cf. 



Lobanov 1998 and references therein). 



We are aware, however, that this offset may be due 
to the positional uncertainties of measurements and lim- 
ited resolution and blending at lower frequency. Future 
high resolution, multi-frequency VLBI observations will 
be critical to clarify this. 

4-3. Core identification 

It is hard to determine the real core in 3C 138 based on 
the apparent source size and spectral index available now. 
From their high resolution polarization sensitive imaging. 



Cotton et al. (1997) have ruled out the possibility that 



Ce was the core and further suggested that Cw harbors 
the central energy source. The non variation in the flux 
density of Ce, discussed in Sect. 3.2, can also be viewed 
as an indication that C^; probably is not related to the 
center of activity. Future observations with better reso- 
lution at lower frequencies (such as the 1.6 GHz space 
VLBI (VSOP) imaging observations (cf. Hirabayashi et al. 
1998 )) and better sensitivity at higher frequencies should 
be able to clearly separate and detect the two compact 
components Cw and Ce and, thus enable us to undoubt- 
edly determine if the spectral index, as suggested by Cot- 
ton et al. (1997), is flatter in component Cw- 

Our new measurements confirm the systematic in- 
creasing in the knot separation with time. But the rate is 
three times slower than the previous one. The detected su- 
perluminal speed, when combined with the observed two- 
sided jet structure in 3C 138, poses very strong constraints 
on the jet kinematics. 

The expected ratio of the flux density between jet and 
counter-jet is given by R = ijz^^)'^^"- Here, intrinsi- 
cally symmetrical two-sided jets with the same speed in 
two opposite directions are assumed. 

The ratio of the flux density of the eastern to west- 



ern lobe seen in the 1.7 GHz image (Cotton et al. 1997) 
is about 16. Using the spectral index between 1.7 and 



5.0 GHz of 0.75 (Fanti et al. 1989), we obtained the prod 
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uct l3cos9 = 0.46, in good agreement with the value esti- 
mated at 5 GHz (Fanti et al. 1989). But this restricts to 



a maximum apparent speed of 1.6c, which is significantly 
smaller than 3.3c. Therefore, the superluminal motion of 
3.3c in the central core region, though much less than 
earlier reported one, is still too high to allow the overall 
counter-jet structure to be seen in terms of the projection 
effect and Doppler boosting of an intrinsically identical 
jet. 

One solution to circumvent this is to assume a non- 
uniform distribution of R for the whole jet since usually 
lobes are observed to move at a much lower speed. This 
can naturally explain the absence of parsec-scale counter- 
jet in the central region. Cotton et al. ( |1997[ ) obtained a 
lower limit to the ratio R of 20 within the central area. Our 
new observations slightly raised this lower limit to 30 at 
5 times RMS noise level. These lower limits are consistent 
with the detected superluminal motion in the inner region 
being due to the relativistic boosting in the jet. 

Such a change in R could be caused by the strong in- 
teraction of the jet with the surrounding medium in the 
immediate region of the core. This is supported by the 
existence of a large amount of ionized gas in the central 
region (Cotton et al. 1997 ). The bending in the jet or a 



helical jet model is a possible explanation, but unlikely 
to affect R significantly in case of 3C 138 (cf. Fanti et 
al. 1989| ). A superluminal speed of 3.3c implies a ratio R 
greater than 360 (using the spectral index of component 
Ce in Table 1), 12 times the lower limit estimated above. 
On the other hand, the possibility that the observed over- 
all asymmetry is intrinsic is still there, but difficult to test 
because of the limited data gathered. The nature of the 
core region is worthy of a more detailed study. 
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